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EOOY  CURRENT  CALCULATIONS  IN  THIN  CONDUCTING 
PLATES  USING  A  FINITE  ELEMENT-STREAM  FUNCTION  CODE* 

K-Y.  Yuan.  J.F.  Abel,  F.C.  .Moon 


Abstract  -  A  stream  function  or  vector  potential  for 
the  current  density  vector  Is  used  to  develop  a  finite 
element  code  for  calculating  Induced  currents  In  thin 
conducting  plates.  While  two-dimensional,  the  code  In- 
eludes  self  field  effects  for  harmonic  fields  for  skin 
depths  on  the  order  of  the  plate  thickness  and  larger. 
The  solution  of  the  resulting  Integro-dl f ferentlal 
equation,  using  a  Galerkln  method,  leads  to  a  complex, 
nonsymmetrlc,  fully  populated  global  matrix.  In 
addition  to  the  current  density,  the  code  also  calcu¬ 
lates  Induced  temperatures  due  to  Joule  heating  and  the 
magnetic  forces  on  the  plate.  The  results  are  compared 
to  Infrared  measurements  of  Induced  currents  in  rec¬ 
tangular  plates.  Extension  of  the  code  to  transient 
problems  using  both  fast  Fourier  transform  and  direct 
integration  methods  Is  In  progress. 

INTRODUCTION 

The  goal  of  ongoing  research  at  Cornell  is  to 
develop  numerical  methods  for  the  coupled  analysis  of 
forces,  currents,  and  stresses  In  thin  elastic 
structures  exposed  to  time-dependent  magnetic  fields. 

In  particular  the  project  Is  directed  toward  time  de¬ 
pendent  magnetic  forces  and  motions  of  elastic  plates 
and  shells.  Soth  applied  mathematical  analysis  and 
finite  element  numerical  computer  techniques  are  being 
used.  Verification  of  the  newly  developed  comouter 
codes  is  carried  out  by  direct  comparison  of  the  calcu¬ 
lated  results  with  experimental  data  and  Is  carried  out 
at  each  stage  of  the  development.  The  experiments 
employ  a  new  Infrared  scanning  technique  to  visualize 
the  induced  eddy  current  patterns  In  the  structures. 

Most  techniques  for  calculating  stresses  due  to 
electromagnetic  fields  require  a  three-dimensional 
treatment  even  If  the  structure  Is  one  or  two- 
dimensional  such  as  a  beam,  plate,  or  shell  [1,2].  In 
this  research  a  stream  function  for  the  Induced  current 
In  the  conducting  structure  has  been  used  which  reduces 
the  magnetics  problem  to  a  two-dimensional  one  and 
allows  one  to  use  conformable  finite  element  meshes  for 
the  magnetic  and  elastic  deformation  fields. 


dimensional  space,  and  constant  across  the  sheet  thick¬ 
ness.  Application  of  theTB1ot-Savart  law  and  the  di¬ 
vergence  theorem  yields  B*  for  the  midplane  of  the 
plate  (z-0): 
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*  has  been  set  to  zero  on  the  boundary  of  the 
plate  In  the  above  derivation.  When  the  plate  is 
multiply  connected,  different  constant  values  of  <s>  will 
be  assigned  on  each  Interior  boundary.  These  boundary 

conditions,  together  with  the  known  3°  and  (3),  are 
then  used  In  (2)  for  the  determination  of  the  stream 
function  >t». 

Formulation 

For  steady  state,  harmonic  current  In  a  flat  plate 
plate,  (2)  may  be  nondlmensional 1  zed  into  the  following 
form 
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in  which  the  coordinates  are  nondlmensional  1  zed  with 

,hB,  iwt 

respect  to  half  the  thickness  (4j0 ;  *  »  • 

8°  .  (.|)ee1“t;  the  magnetic  Reynolds  number  Is 
R  •  wuah2/8*  which  is  related  to  skin  depth  i  through 
a  .  (£.)2;  and  B  is  the  reference  magnetic  field. 

The  finite  element  Galerkin  method  is  used  to 
solve  (4).  a  Is  approximated  globally  and  locally  by 


Flat  plate.  The  derivation  for  the  case  of  a  flat 
plate  has  been  given  by  Moon  [3].  The  details  are  not 
repeated  here;  only  the  essential  relationships  are 
recapitulated. 

If  the  mid-plane  of  the  plate  Is  chosen  as  the  xy 
plane,  one  obtains 

I  «  r  x  (yn)  »  -n  x  7,  (1) 

and 

72*  *  oh  1^-  (8°  *  3*)  .  (2) 

For  the  determination  of  3*.  the  current  stream 
function  >  is  assumed  to  be  defined  throughout  three- 
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In  which  G  Is  the  total  number  of  nodal  points,  E 
denotes  the  Eth  element,  \  are  the  quadratic  global 
Interpol atlori^unctlons  generated  from  the  local  ele¬ 
ment  shape  functions  NE.  Six-node  triangular  elements 

are  used  here.  The  local  element  shape  functions  are 
all  quadratic  In  this  case.  The  element  algebraic 
equations  are 
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Values  of  the  stream  function,  eddy  current, 
temperature  Induced  In -a  half  cycle  of  the  current,  and 
time-averaged  magnetic  pressure  are  produced  as  output. 
The  stream  function  Is  calculated  at  the  nodes  of  the 
finite  element  mesh.  Current,  temperature,  and 
pressure  are  evaluated  at  the  centroid  of  each  element. 
The  stream  function  and  current  are  calculated  in  com¬ 
plex  form.  The  modulus  and  phase  angle  of  the  current 
are  evaluated  in  the  interest  of  spectral  analysis  for 
the  calculation  of  transient  currents. 


The  weighting  function  3j  is  given  by 

c  N^(x.y) 

tf.-U.n)  »  /  - w - j - jtj  dxdy 

J  ,E  C(x-c)2+(y-n)Mr/Z 


(8) 


Because  of  the  numerical  difficulty  associated  with  the 
sharp  variation  of  the  Kernel  function  in  (8),  the 
weighting  function  is  calculated  analytically  within 
the  element  and  numerically  outside  the  element.  (7c) 
is  then  integrated  entirely  by  numerical  quadrature. 

Since  is  quadratic,  six  basic  integrations  with 

nwnerators  1,  x,  y,  x’,  xy,  and  y^  in  (8)  need  to  be 
performed.  To  simplify  the  integration,  each  element 
is  first  .mapped  onto  a  standard  triangle,  which  is 
independent  of  the  coordinates  of  the  nodes  of  the  ele¬ 
ment.  The  integrations  are  then  performed.  The  ex¬ 
pressions  obtained  are  then  used  for  the  calculations 
of  (3)  for  all  the  elements. 

The  resulting  global  matrix  Is  complex,  fully 
populated,  and  nonsymmetrlc  for  general  nonuniform 
meshes. 


CAPABILITIES  of  THE  PROGRAM  ED0Y2 

A  fortran  program  EDDY2  has  been  developed  based 
on  the  formulation  of  (5)-(8).  It  calculates  the  local 
and  nonlocal  solutions  of  stream  function,  eddy  cur¬ 
rent,  temperature,  and  pressure.  As  of  this  writing, 
the  image  solution  and  two-dimensional  graphic  output 
capabilities  remain  to  be  implemented.  Uniform  mag¬ 
netic  field  and  fields  due  to  any  nunber  of  magnetic 
dipoles  can  be  handled.  Magnetic  fields  generated  from 
some  types  of  colls  of  interest  will  be  added. 

Input 

The  geometry  of  the  plate  and  the  description  of 
the  external  magnetic  field  are  the  two  basic  forms  of 
data  needed  by  program  ED0Y2.  The  total  nunbers  of 
nodal  points,  load  cases,  and  elements  need  to  be 
specified.  Coordinate  and  boundary  condition  must  be 
given  for  each  input  node.  Intermediate  nodes  may  be 
generated  for  any  groups  of  nodal  points  that  are  uni¬ 
formly  spaced.  Element  Informaiton  may  also  be  gener¬ 
ated.  Although  only  six-node  triangles  are  included  in 
the  present  version  of  ED0Y2,  the  program  has  been 
structured  so  that  other  types  of  elements  may  be 
added.  Element  group  information  and  the  master  card 
'or  each  element  group  need,  therefore,  to  be  inputed 
too. 


The  program  allows  for  different  orders  of  numeri¬ 
cal  integrations.  Six-  and  seven-point  formula  are  now 
provided.  The  order  may  be  specified  on  the  master 
element  group  card. 


Examples 

Two  problems  analyzed  by  ED0Y2  are  shown  below. 
Figure  1  shows  the  stream  function  contours  for  a  A  to 
1  rectangular  plate  excited  by  a  harmonic  uniform  field 
with  magnetic  Reynolds  number  R  »  0.0012.  Figures  2 
and  2  give  the  stream  function  contours  and  isotherms 
for  a  notched  plate  subjected  to  a  harmonic  uniform 
field  with  R  -  0.001.  Only  half  of  the  plate  is  shown 
In  both  problems  because  of  symmetry.  In  the  analysis, 
however,  the  whole  plates  are  analyzed  without  taking 
advantage  of  the  symmetry  conditions. 

The  CPU  time  is  145  seconds  with  32  elements  in 
the  first  problem  using  the  IBM  370-168  computer.  Most 
of  the  time  Is  used  in  evaluating  the  nonlocal  inte¬ 
gration  terms  in  (6).  The  computational  time  can  obvi¬ 
ously  be  reduced  by  using  the  symmetry  conditions  and 
by  optimizing  the  program.  Also,  the  quadratic  rec¬ 
tangular  element  with  it  much  simpler  experssions  for 
the  nonlocal  integration  terms  will  substantially 
reduce  the  CPU  time. 


Fig.  1  Strecm  function  contours  for 
a  long  rectangular  slate  ex¬ 
cited  b«  a  harmonic  uni  form 
field  (R  =  0.3012). 


REFERENCES 


Ft".  2  Stream  function  contours 
for  a  notched  plate  ex¬ 
cited  by  a  harmonic  uni¬ 
form  field  (ilotchwldth 
■  2h,  S  «  0.001). 


[1]  K.  Miya,  S.  An,  Y.  Ando,  M.  Ohta,  Y,  Suzuki, 
“Application  of  Finite  Element  Method  to  Electro- 
Magneto-Mechanical  Dynamics  of  Superconducting 
Magnet  Coll  and  Vacuum  Vessel,"  Proc,  6th  Symp.  of 
Engineering  Problems  of  Fusion  Research,  No.  I97b,~ 
Publ.  IEEE,  NY  1976,  pp.  927-9$4. 

[2]  E.  B.  Becker,  and  R.  D.  Plllsbury,  “Finite  Element 
Analysis  of  Coupled  Electric,  Magnetic,  and 
Dynamic  Problems,"  Formulations  and  Computational 
Algorltiws  in  Finite  Element  analysis,  K.  J.  Bathe 
et  a T. 7 Ms. HIT  Press', T977',  "ppTTHS9-lQ83. 

£3]  F.  C.  Moon,  "Problems  in  Magneto-Solid  Mechanics,'1 
Mechanics  Today,  S.  Nemat-Nasser,  Ed.,  American 
Academy  of  Mechanics,  Nov.  1977. 


Fig.  3  Isotherms  for  a  notched 
Plate  Excited  by  a  Har¬ 
monic  Uniform  Field 
('!otchv:idth  *  Eh, 

3  *  3.301). 


Report  7,  Part  2 


A  BOUNDARY  INTEGRAL  METHOD  FOR  EDDY  CURRENT  FLOW 
AROUND  CRACKS  IN  THIN  PLATES 


M.A.  Morjaria,  S.  Mukherjee,  and  F.C.  Moon 


presented  at  COMPUMAG-3,  the 

Third  International  Conference  on  the  Computation  of  Magnetic  Fields 

Chicago 

September  1981 


to  be  published  in 
IEEE  Transactions  on  Magnetics 


March  1982 


Unclassified 


StCUKITy  CLASSIFICATION  of  THIS  PAOE  fWhwi  Data  E nlmrtd) 


REPORT  DOCUMENTATION  PAGE 


\ .  REPORT  NUMBER 


7,  Part  2 


2,  GOVT  ACCESSION  NO 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


1.  RECIPIENT'S  catalog  number 


4.  TITLE  (fid  Submit) 


A  BOUNDARY  INTEGRAL  METHOD  FOR  EDDY  CURRENT 
FLOW  AROUND  CRACKS  IN  THIN  PLATES 


5.  TYPE  OF  REPORT  4  PERIOD  COVERED 

Topical  Report 
August  1981 -August  1981 


«.  PERFORMING  ORG.  REPORT  NUMBER 


7.  AuTMORfi; 

M.A.  Morjaria,  S.  Mukherjee  and  F.C.  Moon 


8.  CONTRACT  OR  GRANT  NUMBERf.J 

ONR  Contract  Number 
N00014-79-C-0224 


».  PERFORMING  ORGANIZATION  NAME  ANO  ADDRESS 

Departments  of  Structural  Engineering  and  Theo¬ 
retical  &  Applied  Mechanics,  Cornell  University, 
Ithaca,  NY  14853 


10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  4  WORK  UNIT  NUMBERS 


NR  064-621 


CONTROLLING  OFFICE  NAME  AND  ADORESS 

Director,  Structural  Mechanics  Program,  Material 
Sciences  Division,  Office  of  Naval  Research, 
Arlington,  VA  22217 _ 


12.  REPORT  OATE 

September  1981 


13.  NUMBER  OF  PAGES 
6 


14.  MONITORING  AGENCY  name  4  AODRESS (II  dill, ton t  Itam  Controlling  Ollleo) 


IS.  SECURITY  CLASS,  (o I  thl,  raporlj 


Unclassified 


DECLASSIFICATION/  DOWNGRADING 
SCHEDULE 


i«.  distribution  statement  (ot  thia  Report) 


This  document  has  been  approved  for  public  release  and  sale;  distribution 
unlimited. 


17.  DISTRIBUTION  STATEMENT  (ot  tha  abatract  antarad  In  Block  20,  It  dltlarant  from  Raport) 


IS.  SUPPLEMENTARY  NOTES 


19.  KEY  WOROS  (Continua  on  raaaeaa  atda  it  nacaaamry  and  idantity  by  block  numbar) 

Boundary  element  method,  Eddy  currents,  Nondestructive  testing. 
Numerical  methods.  Plates 


ABSTRACT  (Continua  on  ravaraa  aida  it  nacaaaary  and  Idantity  by  btock  numbar) 

'  A  boundary  element  method  which  employs  a  Green's  function  for  a  crack 
has  been  developed  to  calculate  the  induced  eddy  current  flow  around 
cracks  in  thin  conducting  plates.  The  theoretical  equations  employ  a 
stream  function  for  the  current  density  vector  and  is  equivalent  to  the 
electric  field  vector  potential  method.  A  low  frequency  or  large  skin 
depth  approximation  leads  to  a  Poisson  equation  for  steady  harmonic 
inductor  fields.  Induced  currents  around  a  crack  in  a  square  plate  due  _ 


r?4 ^ 


DD  |  JAN*73  1  47  3ft  EDITION  OF  I  NOV  S3  IS  OBSOLETE 
VN  0102-  LF-  014-  4601 


Unclassified 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Whan  Data  Bntarad) 


b>to  a  uniform  inductor  field  for  various  crack  positions  and  sites  have 
been  calculated*  in  this  papery 

r>The  effect  of  the  relative  position  and  length  of  the  crack,  'J’JP /“acif 
to  the  plate  width,  on  the  eddy  current  density  near  the  tips  of  the  crack 
is  given  special  attention.  These  results  may  be  useful  to  simulate  eo  y 
current  flow  detection  phenomena.^ _ 


S/N  0102-  LF-0H-«60’ 

Unclassified 


SECURITY  CLASSIFICATION  OF  THIS  PAGEC**»n  D.f.  Enter. tf) 


A  BOUNDARY  INTEGRAL  METHOD  FOR  EDDY  CURRENT  FLOW  AROUND  CRACKS  IN  THIN  PLATES 


M.A.  Morjaria,  S.  Mukherjea  and  F.C.  Moon 
Departmanc  of  Theoracical  and  Applied  Mechanics 
Cornell  University,  Ithaca,  Saw  Yorlc 


ABSTRACT 

» 

A  boundary  element  nathod  which  employs  a  Green's  function  for  a  craclc  has  been  developed  to  calculate  the 
induced  addy  current  flow  around  cracks  in  thin  conducting  plates.  The  theoretical  equations  employ  a  stream 
function  for  the  current  density  vector  and  is  equivalent  to  the  electric  field  vector  potential  method.  A  low 
frequency  or  large  skin  depth  approximation  leads  to  a  Poisson  equation  for  steady  harmonic  inductor  fields. 
Induced  currents  around  a  crack  in  a  square  place  due  to  a  uniform  Inductor  field  for  various  crack  positions 
and  slcea  have  been  calculated  in  chis  paper. 

The  effect  of  the  relative  position  and  length  of  Che  crack,  with  respect  to  the  plate  width,  on  the  eddy 
current  density  near  che  tips  of  Che  crack  is  given  special  attention.  These  results  may  be  useful  to  simulate 
eddy  current  flow  detection  phenomena. 


INTRODUCTION 

The  boundary  element  method  (BEM)  (also  called 
che  boundary  integral  equation  method)  has  emerged  as 
an  important  computational  technique  for  eleccrodynamlc 
problems.  Vu  ec  al  [1]  and  Ancelle  ec  al  [2]  have  ad¬ 
dressed  magnetostatic  problems  by  che  BEM  while  Trow¬ 
bridge  [3]  has  considered  problems  by  the  magnetic  po¬ 
tential  method.  Very  recently.  Salon  and  Schneider 
[4J  have  solved  problems  of  eddy  current  flow  in  long 
prismatic  conductors  by  che  BEM  besed  on  the  electric 
potential  approach. 

In  chis  paper,  we  describe  a  powerful  boundary 
element  technique  for  calculating  induced  eddy  current 
flows  in  conducting  places  with  through  cracks  using 
che  electric  potential  approach.  The  BEM  has  the  im¬ 
portant  advantage  that  only  the  boundary  of  a  body 
(.rather  than  che  entire  domain)  needs  to  be  discreclsed 
in  a  numerical  solution  procedure. 

There  have  been  soma  attempts  to  model  eddy  cur¬ 
rent  flow  around  annular  cracks  in  rods  and  in  plates 
by  replacing  cracks  by  slots  (see  for  example  Ref. 

[ 5 ] ) .  However,  we  have  shown  that  the  induced  current 
in  che  vicinity  of  a  crack  leads  to  a  singularity  of 
current  density  at  che  crack  tips  [6,7].  This  high 
concentration  allows  one  to  use  eddy  current  testing 
devices  such  as  active  and  passive  search  coils  to 
decect  che  presence  of  cracks.  It  also  results  In  a 
temperature  hot  spot  which  can  be  detected  by  infrared 
scanning  (6,8).  The  boundary  element  technique  intro¬ 
duced  by  che  authors  [6,7]  and  described  here  allows 
sne  to  model  exactly  the  singular  nature  of  current 


density  at  crack  tips  of  chin  places.  This  technique 
can  handle  any  arbitrary  shape  of  the  place  and  gener¬ 
al  magnetic  fields. 

In  this  paper  we  discuss  application  of  the  BEM 
to  eddy  current  flow  in  a  cracked  square  plate  due  to 
an  uniform  inductor  field  applied  normal  to  the  plate. 
A  number  of  crack  sizes  to  place  size  configurations 
has  been  considered.  Also,  effect  of  the  relative 
position  of  a  crack  tip  to  che  plate  edge  on  the  in¬ 
duced  eddy  current  distribution  has  been  investigated. 

GOVERNING  EQUATIONS 

A  chin  plate  with  a  crack  in  it  is  shown  in  Fig. 
1.  The  plate  is  made  of  a  conducting  material  of 
conductivity  o.  The  plate  boundary  can  be  arbitrary 
and  its  thickness  (uniform)  is  h.  The  thin  line 
crack  is  of  length  2a  and  can  have  arbitrary  ori¬ 
entation  relative  to  the  boundary  of  Che  plate.  The 
coordinate  system  for  the  problem  is  also  shown  in 
Fig.  1.  The  origin  of  coordinates  lies  at  the  center 
of  the  crack  and  at  the  midsurface  of  the  plate. 

An  external,  oscillatory  magnetJc  field,  g°, 
is  applied  which  induces  a  current  density  J  in  the 
plate.  It  is  assumed  that  the  current  density  is 
uniform  across  the  plate  thickness  and  that  the  skin 
depth  (which  Is  inversely  proportional  to  the  square 
root  of  the  frequency)  is  large  compared  to  che  plate 
thickness. 


(7) 


A  scream  function  (or  electric  potential)  formu¬ 
lation  Is  used  in  this  problem.  The  stream  function, 
tj 1*  defined  as 

J  -  7x(>i,k)  -  -k  x  7*  (1) 


This  equation  guarantees  Che  conservation  of  charge 
equation  7>J  »  0  for  charge  free  regions. 

Using  Ohm’s  law  the  governing  differential  equa¬ 
tion  for  che  stream  function  is  obcained  as  [6,7] 


7”i) 


17  <B3+B3> 

»i 


(2) 


In  the  above,  B^  is  the  self  magnetic  field 
due  to  the  current  J.  It  has  been  shown  in  ref. 
[9],  however,  chat  for  a  sinusoidal  applied  field, 
with  che  skin  depch  much  greater  than  the  thickness 
of  the  place,  3^  can  be  neglected  relative  to  the 
applied  field  3  .  This  assumption  simplifies  the 


problem,  and,  with  B¬ 


B3 


(with  i  •  i£T 

■'j  sacis- 


and  u  the  frequency) ,  the  spatial  part  of 
fies  a  two-dimensional  nonhomogeneous  Poisson's 
equation 

72!/  *  tiioBj  •  ffx^Xj) 


(3) 


The  boundary  condition  requires  that  the  currant 
must  be  tangential  to  che  place  boundary.  Thus  tji 
is  required  to  be  constant  on  che  boundaries 
and  jCj .  On  one  boundary,  che  value  of  v  is  set 
to  zero,  while  on  the  other  boundary  j,  «  C  and  C 
is  obtained  from  che  assumption  chat  the  net  flux  flot^* 
ing  through  the  crack  boundary  is  zero.  This  leads  to 
che  condition 


0  J.tds  •  0  (4) 
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where  c  is  an  unit  tangent  to  3C1  and  s  is  the 
distance  measured  along  a  boundary  in  the  anticlock¬ 
wise  sense.  This  formulation  assumes  that  no  current 
flows  across  the  crack  or  crack  tip  and  leads  to  a 
singularity  of  the  J  field  at  a  crack  tip.  This  is 
analogous  to  the  stress  singularity  in  fracture  me¬ 
chanics.  It  is  possible  that  some  leakage  of  current 
occurs  across  a  crack  tip  and  thus  relieves  the  singu¬ 
larity  in  actual  conductors.  Possible  leakage  of 
current  is  not  considered  in  this  paper.  (It  is  noted 
here  that  infrared  scans  of  eddy  current  flow  around 
cracks  do  indeed  show  a  large  Increase  in  temperature 
at  che  crack  tips,  indicating  high  current  density 
at  che  crack  tips  [6].) 

In  suninary,  the  boundary  conditions  on  i,  used 
In  this  formulation,  are 

v  ■  0  on  che  crack  boundary  3C^  (5) 

*  0  on  the  ouC3ide  boundarv  3C.  (6) 
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These  boundary  conditions,  together  with  the 
field  equation  (3) ,  constitute  a  well  posed  oroblem. 


BOUNDARY  ELEMENT  FORMULATION 


Integral  equations 

An  integral  equation  formulacion  for  Poisson’s 
aquation  (3)  can  be  written  as  (Fig.  1)  [6,7] 

2m*>(p)  -  4  k(p,Q)G(0)ds  +  /AK(o,q)f(q)dAa  (8) 
3C, 

This  is  a  single  layer  potential  formulacion 
where  C,  a  source  strength  function  on  the  outside 
boundary,  must  be  determined  from  che  boundary  condi¬ 
tion  on  it  (equation  9) .  The  points  ?  (or  P)  and 
q  (or  Q)  are  source  and  field  points,  respectively, 
with  capital  letters  denoting  points  on  che  boundary 
of  che  body  and  lower  case  letters  denoting  points 
inside  the  body.  The  area  of  the  body  B  is  denoted 
by  A. 

It  has  been  shown  [6]  that  ji  from  equation  (8) 
with  the  following  kernel  satisfies  the  boundary  con¬ 
ditions  (5)  and  (7)  implicidy . 

K(p,q)  -  Re[i(z,?,zo]  (9) 

p(z,7,z  )  -  ln(l-r./E)  -  in(l-r.T)  (10) 

0  / -  1  1 

z  ifiZ 

where  rt  .  _2_ _J 2_  ,  jrj  <,  1 


Re  denotes  the  real  part  of  the  complex  argumene, 
z  zo  are  t*le  sourcl!  and  field  ooinc  coordinates, 

respectively,  in  complex  notation  and  a  suoeroosed 
bar  denotes,  as  usual,  the  comolex  conjugate  of  a  com¬ 
plex  quantity. 

The  remaining  boundary  condition  (6)  on  the  out¬ 
side  surface  is  satisfied  by  using  a  differentiated 
version  of  (8)  and  caking  the  limit  as  ?  inside  3 
approaches  a  point  P  on  30^.  Defining 

H  -  Im(-*±  -  If)  ,  H,  -  -Re(|£  -  %  (U) 

*  jZ  —  34  . 

3z 

the  boundary  condition  (6)  becomes 


where  n^  are  the  components  of  the  unit  outward  norm¬ 
al  to  3C,  at  some  locally  smooth  point  on  it. 


The  current,  J,  at  a  point  inside  the  body  is 
obcained  from  equations  (1)  and  (8). 


and  the  skin  depth 


Discretization  of  equations  and  solution  strategy 

The  outer  boundary  of  the  body,  3C, ,  is  divided 
into  N,  straight  boundary  elements  using  Nb  (tlb  - 
N.,)  boundary  nodes  and  the  interior  of  the  body.  A, 
is  divided  into  n^  triangular  internal  elemencs. 

A  discretized  version  of  equation  (12)  is 

°  ’  lN2-rAS1Hi(PM'C!)ni(PM)G«)d3Q 

+  A^V^V^^q  (13) 

where  is  the  point  P  where  it  coincides  with 

a  node  M  at  a  center  of  a  boundary  segment  on  3C, 
and  £s^  and  £A^  are  boundary  and  internal  elements 
respectively. 

A  simple  numerical  scheme  is  used  in  which  the 
source  strengths  G  are  assumed  to  be  piecewise  uni¬ 
form  on  each  boundary  segmenc  with  their  values  to  be 
determined  at  the  nodes  which  lie  ac  the  centers  of 
each  segment.  Substitution  of  the  piecewise  uniform 
source  strengths  into  equation  (13)  and  carrying  out 
of  the  necessary  integrations,  analytically  and  num¬ 
erically,  leads  to  an  algebraic  system  of  the  type 

(0)  -  (Al(G>  +  (d)  (14) 

The  coefficients  of  the  matrix  [A]  contain 
boundary  integrals  of  the  kernel.  The  vector  (d) 
contains  contributions  from  the  area  integrals  and 
the  vector  (G)  the  unknown  source  strengths  at  the 
boundary  nodes.  The  dimension  of  { G >  depends  only 
on  the  number  of  boundary  elements  on  oCj  and  the 
internal  discretization  is  necessary  only  for  the 
evaluation  of  integrals  with  known  integrands. 

The  solution  strategy  is  as  follows.  The  matrix 
[A]  and  vector  (d)  in  equation  (14)  are  first  eval¬ 
uated  by  using  the  appropriate  expressions  for  the 
kernels  and  the  prescribed  function  f  in  equation 
(3).  Equation  (14)  is  solved  for  the  vector  (G). 

This  value  of  (G)  is  now  used  in  a  discretized  ver¬ 
sion  of  equation  (8)  to  obtain  the  values  of  the 
stream  function  4  at  any  point  p.  Finally,  the 
current  vector  at  any  point  is  obtained  from  equations 
analogous  to  (8) . 


NUMERICAL  RESULTS 


In  the  numerical  computations,  B^  in  Eq.  (13) 
is  assumed  to  be  a  constant.  Eq.  (3)  can  be  non- 
dimensionalized  to  the  form 
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For  the  results  in  this  paper  a  «  2.  A  typical 
mesh  for  the  results  for  example  shown  in  Fig.  2d 
has  48  boundary  segments  uniformly  distributed  along 
the  upper  half  (due  to  symmetry)  of  the  boundary 
of  the  plate.  In  order  to  evaluate  the  known  area 
integral  in  Equation  13,  the  internal  area  quadra¬ 
ture  was  used.  It  took  about  300  c.p.u.  secs  on 
IBM  370/168  to  obtain  the  results  in  Fig.  2d. 

The  equation  (15)  is  identical  to  one  relating 
to  the  torsion  of  shafts.  The  BEM  was  verified  by 
comparing  the  numerical  results  for  the  solution  of 
(15)  in  a  square  plate  without  a  crack  to  known 
analytical  results  for  the  torsion  of  a  shaft.  The 
BEM  method  has  also  been  checked  against  a  finlce 
element  technique  developed  for  eddy  current  prob¬ 
lems  [10]. 

Eddy  current  stream  lines  ( lines)  are  shown 
in  Figs.  2  and  3  for  a  square  plate  with  a  crack  in 
it.  Fig.  2  (a)  -(c)  shows  how  the  stream  lines 
are  affected  by  varying  the  size  of  the  plate  while 
keeping  the  crack  size  same.  Due  to  symmetry  only 
the  upper  half  of  the  plate  is  shown  in  Fig.  2.  Fig. 
2  (d) shows  the  effect  of  moving  the  crack  towards 
one  of  the  plate  edges.  Fig.  3  shows  a  close  up  of 
Che  stream  lines  near  right  crack  tip  for  Fig.  2  (c) . 
The  crowding  of  stream  lines  near  crack  tps  leads 
to  large  gradient  of  4  and  therefore  large  induced 
currents  in  this  region.  The  local  temperature  is 
proportional  to  Che  square  of  the  current  density 
(J.J).  Figure  4  shows  calculated  temperature  scans 
along  a  line  slightly  above  the  crack  (x2  •  .0125) 
for  the  results  shown  in  Fig.  2.  From  Figs.  4  (a)  - 
(c)  one  can  conclude  that  as  the  crack  size  increases 
relative  to  the  plate  size  the  hot  spots  at  crack 
tips  are  more  significant  compared  to  those  at  the 
edges.  The  effect  of  moving  the  crack  near  the  plate 
edge  gives  rise  to  significant  hot  spots  as  shown 
in  Fig.  4  .'d)  and  (c) .  This  becomes  more  apparent 
when  we  look  ac  the  'Eddy  Current  Intensity  Factor' 
defined  below.  It  has  been  shown  [6,7]  that  the 
eddy  current  density  squared  is  inversely  proportion¬ 
al  to  the  distance  r  from  a  crack  tip.  We  can  de¬ 
fine  an  eddy  current  intensity  factor,  39 


where 


(15) 


Table  1  shows  cha  calculated  values  of  M.TT  f or  cha 
cwo  crack  cipa  for  Cha  raaulcs  shown  In  Fig.  2. 

It  la  saan  that  Cha  value  of  raaalns  practical¬ 

ly  conscanc  for  varying  place  sizes.  However  It 
changes  signif lcancly  as  a  crack  cip  la  brought  near 
an  edge  of  cha  place. 


Table  1.  Scress  Incensicy  Faccor  Mill 
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Right 
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Tip 

Left 

Crack 

Tip 

Figures  2,  4 

0.05  1.0 

0.125 

0.125 

(a) 

0.10  1.0 

0.130 

0.130 

(b) 

0.25  1.0 

0.145 

0.145 

(c) 

0.10  0.6 

3.96 

1.30 

(d) 

0.10  0.3 

15.45 

6.93 

(e) 
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